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1.  Introduction 


Explosive  missile  warheads  were  used  by  the  U.S.  Navy  during  the  Persian 
Gulf  War  to  generate  electromagnetic  pulses  (EMP)  for  the  purpose  of 
disabling  electronic  components  such  as  computer  chips.  The  objective  of 
this  action  was  to  neutralize  communications,  command,  control,  and 
intelligence  (C3I)  systems.  Because  of  the  Navy’s  air  attack  role,  particular 
emphasis  was  given  to  disabling  air  defense  electronic  systems. 

In  this  report,  equations  describing  the  attenuation  of  long-wavelength 
electromagnetic  radiation  by  water  fogs  and  rain  are  briefly  reviewed  and 
explained;  the  assumption  is  made  that  aerosol  droplets  are  spherical  and  that 
irradiance  of  the  beam  is  in  the  linear  regime  (i.  e.,  the  irradiance  is  not  high 
enough  to  produce  nonlinear  effects  such  as  gas  breakdown  and  explosive 
vaporization  of  the  droplets). 

Calculations  of  water  fog  attenuation  coefficients  are  made  using  the 
Rayleigh  small  particle  approximation.  An  error  analysis  of  this  approximation 
is  made  by  comparing  Rayleigh  approximation  calculations  of  absorption 
efficiency  with  exact  Mie  theory  calculations  of  extinction  efficiency.  The 
Rayleigh  approximation  is  appropriate  for  fog  with  drop  diameters  up  to 
100  pm  for  frequencies  less  than  50  GHz.  Comparisons  are  made  of  Rayleigh 
absorption  efficiencies  and  Mie  extinction  efficiencies  for  rain  with  water 
drops  up  to  5.0  mm;  for  frequencies  above  1  GHz,  the  use  of  the  Rayleigh 
approximation  is  not  appropriate.  For  various  rain  rates  and  temperatures, 
numerical  integration  of  the  Mie  extinction  efficiency  is  used  to  compute  the 
attenuation  coefficient  for  rain  with  a  Marshall-Palmer  (MP)  drop  size 
distribution.  These  results  are  compared  with  those  given  by  a  power  law 
relationship  with  coefficients  given  by  Olsen,  Rogers,  and  Hodge  (ORH);  this 
relationship  was  used  in  the  Near  Millimeter  Wave  (NMMW)  Module  of  the 
Electro-Optical  Systems  Atmospheric  Effects  Library  (EOSAEL).  The  ORH 
relationship  gives  very  good  results  at  temperatures  where  the  coefficients 
were  determined  (from  -20  to  20  °C).  As  temperature  increases  above  20  °C, 
the  exact  Mie  calculations  deviate  increasingly  from  the  results  of  ORH.  The 
use  of  the  ORH  equation  is  not  recommended  above  20  °C.  The  exact  Mie 
computations  are  recommended;  this  approach  can  be  performed  at  any 
temperature  and  can  be  quickly  done  on  a  personal  computer. 


5 


2.  Defining  and  Computing  the  Attenuation  Coefficients 
With  Mie  Theory 

D.R.  Brown  [1]  provided  equations  to  help  compute  the  attenuation 
coefficient  for  clouds,  fog,  ice  fog,  and  rain.  A  discussion  of  issues  related  to 
these  natural  obscurants  was  given.  A  more  detailed  treatment  of  the  theory 
is  given  here  but  only  spherical  particles  will  be  considered. 

The  number  density  (i.  e.,  number  per  unit  volume)  of  spherical  particles  with 
diameters  between  D,  and  D2  is  given  by  integration  of  the  size  distribution 
function  N  with  respect  to  the  sphere  diameter  D  in  the  following  equation: 

An  (D„D2)  =  j  N  dD.  (1) 

D, 

If  Dj  =  0  and  D2  =  °°,  the  total  number  density  is  obtained.  The  extinction 
cross  section  [2]  aext  is  usually  a  function  of  D  and  the  attenuation  coefficient 
aext  is  as  follows: 


=  J  Noext  dD. 
0 


(2) 


In  the  absence  of  multiple  scattering,  the  significance  of  the  attenuation 
coefficient  aext  is  the  irradiance  I  of  a  wave  propagating  in  the  z  direction  is 
given  as  a  function  of  z  by  the  following  expression: 

I(z)  =  I(z  =  0)  e<_a“'z),  (3) 


where  the  choice  of  z  =  0  is  arbitrary.  The  attenuation  coefficient  aext  has  the 
unit  of  (length)-1.  In  engineering  literature,  an  attenuation  coefficient  a  with 
the  units  of  decibels  (dB)  per  unit  length  (usually  km)  is  commonly  used.  The 
way  to  express  the  relation  between  a  and  aext  which  holds  for  any  unit  of 
length  is  as  follows: 


a  = 


10  dB 
ln(10) 


a 


ext 


»  [4-343  dB]  a„. 


(4) 
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The  reason  for  introducing  a  is  not  always  clear.  However,  recalling  that 
exp(-x)  =  10'x/ln(10)  allows  equation  (3)  to  be  rewritten  with  the  help  of 
equation  (4)  as  follows: 


I  (z) 

I  (z  =  0) 


-az 

JQlOdB 


(5) 


Equation  (5)  is  a  simple  way  to  determine  how  many  orders  of  magnitude  an 
electromagnetic  wave  is  attenuated  at  a  particular  range  z;  this  is  the  reason 
for  introducing  a.  Note  that  10  dB  =  1  bel  (B);  the  convention  is  to  use  dB 
rather  than  B. 

The  power  attenuated  by  a  particle,  Pext,  is  equal  to  the  sum  of  the  power 
scattered  by  the  particle  and  the  power  absorbed  by  the  particle.  Absorbed 
power  refers  to  the  heat  generated  in  the  particle  due  to  Joule  heating  under 
irradiation  by  an  electromagnetic  wave. 

The  extinction  cross  section  is  defined  as  follows: 


(6) 

1inc 

where  Iinc  is  the  irradiance  of  the  incident  beam.  The  extinction  efficiency, 
Qext,  is  defined  as  follows: 


(7) 


where  Pinc  is  the  power  of  the  incident  beam  passing  through  the  area  of  the 
particle  if  it  is  not  there.  Since  the  following  is  true: 


then 


^inc  ^inc 


r  d2^ 

7t  — 


°ex«  =  Q 


ext 


( 

n — 


(8) 


(9) 
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and  substituting  into  equation  (2)  gives  the  following: 


=  7  j  CL,  d2n  dD.  (10) 

4  o 

A  small  particle  approximation  for  Qext  is  used  for  fog,  allowing  equation  (10) 
to  be  integrated  analytically  without  specifying  a  particular  particle  size 
distribution.  For  rain,  equation  (10)  is  evaluated  by  numerically  integrating 
over  a  particular  particle  size  distribution;  Qext  is  computed  by  using  the 
following  Mie  theory  expression:  [2] 

Q„,  =  2r  £  (2n  +  l)  Re(a„+b„).  (11) 

A  N  =  1 

The  an  and  bn  are  referred  to  as  scattering  coefficients  [2]  and  X  is  the  Mie 
size  parameter  which  is  defined  by  the  following  relation: 

X  =  Jt|,  (12) 

where  X  is  the  wavelength  of  the  incident  electromagnetic  wave  in  the 
surrounding  medium.  Before  computing  Qext  using  either  Mie  theory  or  a 

A 

small  particle  approximation,  it  is  necessary  to  obtain  N  the  complex 
refractive  index  of  the  sphere’s  medium  relative  to  the  surrounding  medium; 
this  is  a  function  of  A..  For  a  harmonic  wave  with  time  dependence  of 
exp(-icot),  N  is  written  as  follows: 

N  =  Re(N)  +  ilm(N),  (13) 

A  A. 

where  Re(N)  and  Im(N)  are  both  greater  than  zero.  In  the  appendix,  a 

A 

method  for  computing  N  for  water  in  the  frequency  range  1  -  1000  GHz  is 
discussed;  this  method  is  used  for  both  rain  and  fog. 
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3.  Small  Particle  (i.  e.,  Rayleigh)  Approximation  for  otext 


In  a  particle  size  range  where  a  simple  analytical  approximation  for  the 
extinction  efficiency  Qext  is  known,  equation  (10)  can  be  analytically 
integrated  to  produce  an  expression  for  the  attenuation  coefficient.  Bohren 
and  Huffman  [2]  discuss  approximations  which  hold  if  spherical  particles  are 
small,  producing  the  following  result: 

Qabs  =  4XIm(Z),  (14) 


where 


A  _ 

N2  -  1 


/v  0  • 

N2  +  2 


(15) 


According  to  Bohren  and  Huffman,  [2]  equation  (14)  holds  if  the  following  is 
true: 

|x3Im(Z)  «  1.  (16) 

If  equation  (14)  is  a  good  approximation  and  if  the  scattering  efficiency  is 
much  less  than  the  absorption  efficiency,  then  the  following  is  true: 

Qext  -  4XIm(Z).  (17) 

Substituting  equation  (17)  into  equation  (10)  produces  the  following: 

aext  =  y-  Im  (Z)  N  dD.  (18) 

A  0  6 

The  mass  of  water  Mw  in  a  given  volume  of  water  Vw  is  as  follows: 


Mw  —  Pw^w’ 

where  pw  is  the  density  of  the  water. 


(19) 
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The  following  relation  is  also  useful: 


Mw  =  CVa,  (20) 

where  C  is  the  water  content  (i.  e.,  water  concentration  =  mass  of  water  per 
unit  volume  of  air)  and  Va  is  the  volume  of  air. 

Setting  equation  (19)  equal  to  equation  (20)  and  rewriting  the  result  produces 
the  following: 


Vw  =  _C 
Va  Pw 


and  since  the  following  is  also  true: 


N  dD, 


equation  (18)  can  be  written  as  follows: 


htt  t  /  ry  \  0 

a~.  =  —  Im(Z)  — 
X  p 


ext 


W 


(21) 


(22) 


(23) 


Since  the  index  of  refraction  of  the  air  surrounding  the  water  droplets  does 
not  differ  greatly  from  vacuum,  a  reasonable  approximation  is  that  X  —  C/f, 
where  C  =  2.9979xl010  cm/s  is  the  speed  of  light  in  vacuum  and  f  is  the 
frequency.  This  approximation  allows  equation  (23)  to  be  rewritten  as 
follows: 


67if  T  C 

«ext  =  -TT  Im(z)  — •  (24) 

Pw 

In  the  Rayleigh  approximation,  aext  is  proportional  to  C,  the  concentration  of 
water  in  the  atmosphere.  Initial  conclusions  indicate  that  aext  is  also 
proportional  to  the  frequency  f;  however,  Z  is  a  function  of  frequency,  so  aext 
is  not  exactly  proportional  to  f. 


4.  Computer  Calculations  of  a  for  Fog  Using  the  Rayleigh 
Approximation 

Liebe  et  al.  [3]  have  given  Rayleigh  calculations  of  a  for  water  for 
f  =  0  -  1000  GHz  at  T  =  20  °C  with  C  =  1  g/m3.  Liebe's  published  results  can  be 
used  as  a  benchmark  to  verify  the  accuracy  of  attenuation  calculations 
performed  in  this  investigation. 

For  computer  calculations,  it  is  necessary  to  rewrite  equations  (23)  and  (24)  in 
nondimensional  form.  If  aext  has  units  of  (km)-1,  X  has  units  of  cm,  f  has  units 
of  GHz,  C  has  units  of  g/m3,  and  pw  has  units  of  g/cm3,  then  the 
nondimensional  expressions  of  equations  (23)  and  (24)  appropriate  for 
computer  programming  or  use  with  an  electronic  calculator  are  as  follows: 


0.6ll  I  C' 

v  Im(z)  n> 

A  pw 

(25) 

0.6rcf'  C 

-  Im(Z)  — 

29.979  o' 

(26) 

where 

a'  =  acxt/(km)-> 

X'  =  X/cm 
f'  =  f7GHz 
C  =  C/[g/m3] 

Pw  =  pw/[g/cm3]. 

If  a  has  the  units  of  dB/km,  equation  (4)  may  be  rewritten  as  follows: 


a'  = 


10 


ln(10) 


a. 


(27) 


where  a'  =  oc/[dB/km], 

A  computer  code  based  on  equations  (26)  and  (27)  was  written  to  see  if  the 
results  presented  in  the  work  of  Liebe  et  al.  [3]  for  a  as  a  function  of  f  could 
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be  reproduced.  It  was  necessary  to  compute  the  complex  index  of  refraction 
of  water  as  a  function  of  frequency.  This  was  done  using  the  double-Debye 
formulation  given  by  Liebe  et  al.  [3]  A  review  of  the  double-Debye 
formulation  is  given  in  appendix  A. 

The  results  of  these  calculations  of  a  as  a  function  of  f  are  compared  in 
table  1  with  those  given  by  Liebe  et  al.  [3]  Pruppacher  and  Klett  [4]  state  that 
the  water  content  of  most  fogs  is  less  than  0.1  g/m3.  Therefore,  a  fog  with 
C  =  1  g/m3  is  a  very  heavy  fog;  values  of  a  in  most  fogs  would  be  an  order  of 
magnitude  smaller  than  reported  for  C  =  1  g/m3.  A  water  content  of  C  =  1  g/m3 
represents  a  worst  case,  a  is  also  given  for  C  =  0.1  g/m3,  a  more  typical  case. 


Table  1.  Attenuation  coefficient  a  as  a  function  of  frequency  in  the  small 
particle  approximation  with  pw  =  1  g/cm3 _ 


f/GHz 

Liebe  etal. 

For  C  =  1  g/m3 
a/[dB/km] 

This  work 

For  C  =  1  g/m3 
a/[dB/km] 

This  work 

For  C  =  0. 1  g/m3 
a/[dB/km] 

300 

15.0 

15.801 

1.5801 

400 

20.3 

20.641 

2.0641 

500 

25 

25.197 

2.5197 

800 

38 

36.830 

3.6830 

1000 

43 

42.598 

4.2598 

Liebe  et  al.  [3]  also  shows  curves  of  a  as  a  function  of  frequency  for 
f  =  20  -  1000  GHz  for  temperatures  in  a  range  from  -20  to  +20  °C.  These  curves 
are  reproduced  in  figure  1.  Figure  2  indicates  the  same  curves  for 
f  =  1  -  350  GHz,  showing  more  detail  in  the  lower  frequencies.  At  frequencies 
greater  than  200  GHz,  attenuation  caused  by  water  fog  and  clouds  increases 
with  temperature.  At  frequencies  less  than  50  Hz,  the  attenuation  decreases 
with  temperature  (i.  e.,  the  effect  of  uniformly  heating  water  fog  and  clouds  at 
frequencies  less  than  50  GHz  increases  the  propagation).  Figure  3  shows 
detailed  results  for  a  at  microwave  frequencies  0.1-8  GHz. 


Figure  1.  Attenuation  coefficient  of  water  fog  versus  frequency  for 
f  =  20  —  1000  GHz;  computed  using  the  Rayleigh  approximation  for  water 
concentration  of  1  g/m3  with  T  =  -20  to  20  °C. 
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ATTENUATION  COEF.  ALPHA  /.  (dB/km) 


Figure  2.  Attenuation  coefficient  of  water  fog  versus  frequency  for 
f  =  1  —  350  GHz;  computed  using  the  Rayleigh  approximation  for  water 
concentration  of  1  g/m3  with  T=  -20  to  20  °C. 
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5.  Sample  Calculations  of  Intensity  Attenuation  Using  the 
Attenuation  Coefficient  Curves 

Once  the  attenuation  coefficient  a  is  determined  for  a  fog,  cloud,  or  rain, 
equation  (5)  can  be  used  to  determine  the  attenuation  of  intensity  at  a  given 
range. 

For  example,  if  T  =  20  °C  and  f =  400  GHz,  figure  1  shows  that  a  heavy  fog  with 
C  =  1  g/m3  has  an  attenuation  coefficient  a  approximately  equal  to  20  dB/km. 
Since  a  is  proportional  to  the  water  concentration  in  the  fog,  C  =  0.8  g/m3;  the 
resulting  value  of  a  is  [(20  dB/km)/(l  g/m3)](0.8  g/m3)  (i.  e.,  a  =  16  dB/km). 

At  a  range  of  5  km,  az  =  80  dB;  equation  (5)  then  gives  the  result  that 
I(z)/I(z=0)  =  10'8  so  that  the  electromagnetic  wave  would  be  attenuated  by 
eight  orders  of  magnitude,  a  significant  amount.  For  a  more  typical  fog  under 
the  same  conditions  with  C  =  0.1  g/m3,  az  =  10  dB  and  I(z)/I(z=0)  =  10_1;  the 
electromagnetic  wave  would  be  attenuated  by  only  one  order  of  magnitude. 
The  attenuation  of  a  fog  is  extremely  dependent  on  range  and  water  content. 

The  attenuation  of  microwaves  with  frequencies  less  than  10  GHz  is  much 
less  than  for  higher  frequencies.  If  T  =  20  °C,  C  =  1  g/m3,  and  f  =  8  GHz, 
figure  3  shows  that  a  is  approximately, equal  to  0.04  dB/km.  If  C  =  0.8  g/m3, 
then  a  =  0.032  dB/km;  and  at  a  range  of  5  km,  az  =  0.16  dB.  Equation  (5)  then 
gives  the  result  that  I(z)/I(z=0)  =  10-°-016  =  0.964.  The  beam  is  attenuated  by 
less  than  4  percent  compared  to  eight  orders  of  magnitude  for  a  beam  with 
f  =  400  GHz. 

For  the  more  typical  fog  with  C  =  0.1  g/m3,  az  =  0.020  dB,  and 
I(z)/I(z=0)  =  10-°-0020  =  0.995;  the  attenuation  is  less  than  1  percent  at  a  range 
of  5  km.  At  shorter  ranges,  the  attenuation  is  negligible. 
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Figure  3.  Attenuation  coefficient  of  water  fog  versus  frequency  for 
f  =  0.1  —  8  GHz;  computed  using  the  Rayleigh  approximation  for  water 
concentration  of  1  g/m3  with  T  =  -20  to  20  °C. 
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6.  Analysis  of  Rayleigh  Attenuation  Coefficient  Calculations 


The  Rayleigh  approximation  used  for  these  calculations  is  of  great  interest 
because  it  allows  the  attenuation  coefficient  to  be  calculated  for  fog  or  a 
cloud  without  specifying  a  particular  size  distribution.  The  Rayleigh  results 
are  independent  of  the  particular  size  distribution.  The  only  requirement  is 
that  the  particles  be  small. 


Accuracy  of  Rayleigh  attenuation  coefficient  calculations  can  be  determined 
for  a  particular  size  distribution  by  numerically  integrating  equation  (10);  a 
more  general  qualitative  approach  is  to  compare  the  Rayleigh  absorption 
efficiency  with  the  Mie  extinction  efficiency  over  a  size  range  for  a  given 
temperature  and  frequency.  This  is  done  in  the  next  series  of  figures  for  water 
droplets  with  diameters  of  0.5  -  100  pm  at  temperatures  of  -20  and  20  °C. 
Figures  4  through  13  are  for  frequencies  of  10,  50,  300,  600,  and  1000  GHz. 
Before  analyzing  the  results,  equation  (15)  is  rewritten  as  follows: 


X  * 


~  (1.04793)(1<T5) 


(28) 


Equation  (28)  shows  that  an  increase  in  frequency  corresponds  to  an  increase 
in  the  Mie  particle  size  parameter,  and  X  increases  to  about  1  when 
f  =  1000  GHz  for  D  =  100  pm. 

Figures  4  and  5  show  close  agreement  between  the  Rayleigh  and  Mie  results 
at  10  GHz  with  more  extinction  predicted  for  the  lower  temperature.  Figures  6 
and  7  indicate  essentially  the  same  results  for  50  GHz,  but  figures  8  through 
13  show  that  the  Rayleigh  results  increase  in  error  as  the  frequency  increases 
up  to  1000  GHz,  at  which  point  X  ~  1  when  D  =  100  pm.  However,  the  error  at 
1000  GHz  in  this  size  range  is  still  less  than  an  order  of  magnitude.  For  fog 
with  maximum  diameters  of  100  pm,  Rayleigh  theory  is  exact  at  microwave 
frequencies  less  than  10  GHz,  within  a  few  percent  at  50  GHz,  and  gives  better 
than  order  of  magnitude  predictions  up  to  1000  GHz.  There  is  increasing  error 
with  the  Rayleigh  approximation  as  frequency  and  particle  size  increase. 

The  Mie  extinction  efficiency  decreases  with  increasing  temperature  for  f  =  10 
and  50  GHz  but  increases  with  temperature  for  f  =  300,  600,  and  1000  GHz. 
This  agrees  with  the  Rayleigh  result  that  attenuation  decreases  with 
increasing  temperature  for  microwaves  propagating  through  water  fog. 
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Figure  4a.  Raleigh  absorption  efficiency  and  Mie  extinction  efficiency  versus  water 
drop  diameter  for  fog  with  D  =  0.5  —  100  #tm,  f  =  10  GHz,  and  T  =  -20  °C. 
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Figure  4b.  Rayleigh  absorption  efficiency  and  Mie  extinction  efficiency  versus  water 
drop  diameter  for  fog  with  D  =  0.5  —  100  pm,  f  =  10  GHz,  and  T  =  20  °C. 
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Figure  5a.  Raleigh  absorption  efficiency  and  Mie  extinction  efficiency  versus  water 
drop  diameter  for  fog  with  D  =  0.5  —  100  jon,  f  =  50  GHz,  and  T  =  -20  °C. 
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Figure  5b.  Rayleigh  absorption  efficiency  and  Mie  extinction  efficiency  versus  water 
drop  radius  for  fog  with  R  =  0.5  —  50  f  =  50  GHz,  and  T  =  20  °C. 
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Figure  6b.  Rayleigh  absorption  efficiency  and  Mie  extinction  efficiency  versus  water 
drop  diameter  for  fog  with  D  =  0.5  —  100  #im,  f  =  300  GHz,  and  T  =  20  °C. 
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Figure  7a.  Rayleigh  absorption  efficiency  and  Mie  extinction  efficiency  versus  water 
drop  diameter  for  fog  with  D  =  0.5  —  100  fim,  f  —  600  GHz,  and  T  =  -20  C. 
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Figure  7b.  Rayleigh  absorption  efficiency  and  Mie  extinction  efficiency  versus  water 
drop  diameter  for  fog  with  D  =  0.5  —  100  jun,  f  =  600  GHz,  and  T  =  20  °C. 
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Figure  8a.  Rayleigh  absorption  efficiency  and  Mie  extinction  efficiency  versus 
water  drop  diameter  for  fog  with  D  =  0.5  —  100  pin,  f  =  1000  GHz,  and 
T  =  -20  °C. 
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Figure  8b.  Rayleigh  absorption  efficiency  and  Mie  extinction  efficiency  versus 
water  drop  diameter  for  fog  with  D  =  0.5  —  100  fim,  f  =  1000  GHz,  and 
T  =  20  °C. 
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7.  Attenuation  Caused  by  Rain  With  an  MP  Size 
Distribution 

Section  6  demonstrated  that  for  fog  with  maximum  diameters  of  100  |im, 
Rayleigh  theory  is  essentially  exact  at  microwave  frequencies  of  1  -  10  GHz. 
This  was  done  by  comparing  the  Rayleigh  absorption  efficiency  with  the  Mie 
extinction  efficiency  over  a  size  range  of  spherical  water  droplets  up  to  100  pm 
in  diameter.  This  section  will  discuss  the  sphericity  of  rain  drops;  in  addition, 
it  will  be  determined  whether  Rayleigh  approximation  can  be  used  for  rain 
drops  in  the  microwave  frequency  range  of  1  -  10  GHz. 

Evidence  suggests  that  rain  drops  can  be  significantly  nonspherical. 
Therefore,  it  is  common  to  introduce  an  equivalent  diameter,  D0,  defined  as 
the  diameter  of  a  sphere  with  the  same  volume  as  the  nonspherical  drop.  This 
approach  was  taken  by  Pruppacher  and  Klett  who  gave  a  succinct  discussion 
of  rain  drop  shapes: 

When  falling  at  terminal  velocity,  drops  are  nearly  perfect  spheres  if 
D0  <  280  pm.  Larger  drops  are  slightly  deformed  and  resemble  oblate 
spheroids  if  280  pm  <  D0  <  1000  pm.  For  D0  >  1000  pm,  the 
deformation  becomes  large  and  the  drops  resemble  oblate  spheroids 
with  flat  bases.  Drops  larger  than  about  10  mm  in  diameter  are 
hydrodynamically  unstable  and  break  up,  even  in  a  laminar  air 
stream.  [4] 

Pruppacher  and  Klett’s  discussion  shows  that  a  nearly  exact  description  of 
rain  drops  is  likely  to  be  extremely  difficult,  requiring  full  time  research  by 
specialists  who  study  the  production  and  mechanics  of  nonspherical  rain 
drops.  Appropriate  electromagnetic  scattering  algorithms  which  can  compute 
Qext  for  nonspherical  drops  are  also  required.  The  results  are  orientation 
dependent  and  require  averaging  over  angles  which  describe  the  orientation. 
The  situation  is  even  worse  if  the  drop  orientations  are  not  isotropic  but  have 
preferred  orientations.  With  all  these  considerations  and  difficulties  in  mind, 
it  becomes  obvious  why  the  assumption  of  a  simple  size  distribution  of 
equivalent  spheres  is  usually  made.  Therefore,  the  given  analysis  of 
attenuation  by  rain  will  assume  spherical  drops. 

The  question  is  whether  the  Rayleigh  approximation  can  be  used  for  rain 
drops  by  comparing  the  Rayleigh  absorption  efficiency  with  the  Mie 
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extinction  efficiency  over  a  size  range  of  spherical  rain  drops  up  to  2.5  mm  in 
the  microwave  frequency  range  of  1  -  10  GHz.  This  is  done  in  figures  14 
through  18  for  frequencies  of  1,  3,  5,  and  10  GHz.  These  figures  show  that  if 
f  =  1  GHz,  the  Rayleigh  approximation  is  within  20  percent;  if  f  =  3  GHz,  the 
Rayleigh  result  is  still  within  an  order  of  magnitude.  This  is  not  true  as 
frequencies  increase  past  5  GHz.  Failure  to  achieve  order  of  magnitude 
agreement  over  the  range  of  interest  resulted  in  the  investigation  of  the 
numerical  integration  of  equation  (10): 

=  7  ]  Q„,D2N  dD.  (10) 

4  0 

D  represents  an  equivalent  sphere  diameter  if  the  drops  are  nonspherical.  If  a 
size  distribution  of  the  rain  drops  is  known  or  assumed,  Qext  can  be  calculated 
using  Mie  theory  and  equation  (10)  can  be  numerically  integrated. 

Oguchi  [5]  discusses  of  the  issues  related  to  nonspherical  drops.  Discussion 
is  made  of  two  commonly  used  size  distributions  of  equivalent  spheres: 
1)  Laws  and  Parsons  and  2)  MP.  Oguchi  performed  the  numerical  integration 
of  these  two  distributions  and  thereby  produced  curves  of  specific 
attenuation  (i.  e.,  dB/km)  as  a  function  of  frequency  with  f  =  1  —  1000  GHz  at 
T  =  20  °C.  Curves  were  produced  for  rain  rates  R  =  2.5,  12.5,  50,  and  150  mm/h. 
According  to  the  National  Weather  Service  Weather  Radar  Code,  [6]  rain 
rates  are  classified  as  light  (0  -  3.8  mm/h),  moderate  (3.8  -  20  mm/h), 
strong  (20  -  40  mm/h),  very  strong  (40  -  76  mm/h),  intense  (76  -  152  mm/h), 
and  extreme  (greater  than  152  mm/h). 
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Figure  9.  Rayleigh  absorption  efficiency  and  Mie  extinction  efficiency  versus  water 
drop  diameter  for  rain  with  D  =  0.05  —  5.0  mm,  f  =  1  GHz,  and  T  =  20  °C. 
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Figure  10.  Rayleigh  absorption  efficiency  and  Mie  extinction  efficiency  versus  water 
drop  diameter  for  rain  with  D  =  0.05  —  5.0  mm,  f  =  3  GHz,  and  T  =  20  °C. 
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Figure  11.  Rayleigh  absorption  efficiency  and  Mie  extinction  efficiency  versus  water 
drop  diameter  for  rain  with  D  =  0.25  —  5.0  mm,  f  =  5  GHz,  and  T  =  20  °C. 
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Figure  12.  Rayleigh  absorption  efficiency  and  Mie  extinction  efficiency  versus  water 
drop  diameter  for  rain  with  D  =  0.05  —  5.0  mm,  f  =  5  GHz,  and  T  =  20  °C 
(scaled  differently  to  facilitate  comparison  with  figure  13). 


36 


EFFICIENCY 


R  /  mm 


Figure  13.  Rayleigh  absorption  efficiency  and  Mie  extinction  efficiency  versus  water 
drop  diameter  for  rain  with  R  =  0.025  —  2.5  mm,  f  =  10  GHz  and  T  =  20  °C. 
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At  Army  Research  Laboratory  (ARL),  a  code  employing  numerical  integration 
of  equation  (10)  has  been  developed;  Oguchi’s  specific  attenuation 
curves  [5]  for  the  MP  size  distribution  are  reproduced  in  figure  19. 
Comparison  of  Oguchi’s  original  and  the  reproduction  shows  close 
agreement.  _ 

Brown  [1]  makes  the  comment  that  full  Mie  scattering  calculations  are  too 
cumbersome  for  inclusion  in  the  EOSAEL.  As  an  alternative  procedure, 
Brown  uses  the  following  power  law  relationship  developed  by  ORH:  [7] 

a  =  aR\  (29) 

where  a  and  b  were  determined  for  various  size  distributions,  including  the 
MP  distribution  used  in  the  numerical  calculations  of  this  report.  The  ORH 
coefficients  (a  and  b)  are  listed  in  tabular  form  as  a  function  of  frequency  for 
41  frequencies  between  1  and  1000  GHz  at  temperatures  of -10,  0,  and  20  C. 
The  ORH  MP  coefficients  for  T  =  20  °C  were  used  to  compute  ORH 
attenuation  coefficients  for  comparison  with  the  integrated  Mie  theory  results 
of  figure  19.  The  results  of  this  comparison  are  shown  in  figure  20.  The  ORH 
results  shown  at  the  41  discrete  frequencies  of  their  table  agree  with  the 
results  given  in  figure  19  for  rain  rates  less  than  or  equal  to  50  mm/h.  More 
deviation  is  seen  for  the  highest  rain  rate  of  150  mm/h.  Figure  21  shows  the 
same  comparison  except  that  only  frequencies  less  than  10  GHz  are 
considered  and  the  same  trend  is  observed  (i.e.,  better  agreement  is  obtained 
for  the  lower  rain  rates).  These  computations  are  repeated  in  figures  22 
through  24  for  T  -  0  °C  with  essentially  the  same  result  (i.e.,  better  agreement 
at  the  lower  rain  rates).  In  EOSAEL  82  to  84,  Brown  [1]  used  the  ORH  MP 
coefficients  [7]  for  a  rain  rate  of  4  mm/h,  which  proved  to  be  appropriate.  The 
main  problem  with  depending  on  the  ORH  approach  is  that  the  coefficients 
are  only  available  for  T  =  -10,  0,  and  20  °C;  however,  the  full  Mie  scattering 
calculations  are  not  cumbersome  and  are  easy  to  produce. 
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Figure  14a.  Attenuation  coefficient  a  of  rain  versus  frequency  for 
f  =  1  —  1000  GHz;  computed  for  rain  rates  R  =  2.5  —  150  mm/h  using  numerical 
integration  of  the  MP  distribution  with  T  =  20  °C. 
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ATTENUATION  COEF.  a  /  (dB/km) 


Figure  14b.  Attenuation  coefficient  a  of  rain  versus  frequency  for 
f  =  1  —  1000  GHz;  compares  numerical  integration  of  the  MP  distribution 
(T  =  20  °C)  with  discrete  points  given  by  the  ORH  equation  (T  =  20  °C). 
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Figure  14c.  Attenuation  coefficient  a  of  rain  versus  frequency  for 
f  =  0.1  —  10  GHz;  compares  numerical  integration  of  the  MP  distribution 
(T  =  20  °C)  with  discrete  points  given  by  the  ORH  equation  (T  =  20  °C). 
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Figure  15a.  Attenuation  coefficient  a  of  rain  versus  frequency  for 
f  =  1  _  1000  GHz;  computed  for  rain  rates  R  =  2.5  —  150  mm/h  using  numerical 
integration  of  the  MP  distribution  with  T  =  0  °C. 
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Figure  15b.  Attenuation  coefficient  a  of  rain  versus  frequency  for 
f  =  1  —  1000  GHz;  compares  numerical  integration  of  the  MP  distribution 
(T  =  0  °C)  with  discrete  points  given  by  the  ORH  equation  (T  =  0  °C). 
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ATTENUATION  COEF.  a  /  (dB/km) 


Figure  15c.  Attenuation  coefficient  a  of  rain  versus  frequency  for 
f  =  0.1  —  10  GHz;  compares  numerical  integration  of  the  MP  distribution 
(T  =  0  °C)  with  discrete  points  given  by  the  ORH  equation  (T  =  0  °C). 
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Deviation  with  increasing  temperature  is  demonstrated  by  figures  25  and  26; 
these  show  integrated  Mie  theory  attenuation  coefficients  for  T  =  40  and 
60  °C.  In  these  figures,  the  ORH  attenuation  coefficients  at  20  °C  are  shown 
for  comparison.  For  the  lower  rain  rates,  figure  25  shows  the  Mie  attenuation 
at  40  °C  becomes  lower  than  that  given  by  ORH  at  20  °C.  However,  agreement 
for  the  higher  rain  rates  is  better. 

Use  of  ORH  coefficients  is  a  good  method  for  avoiding  numerical  integration. 
However,  once  a  code  is  obtained,  it  is  quick  and  easy  to  find  attenuation 
coefficients  by  performing  numerical  integration  of  the  Mie  theory  over  the 
size  distribution  of  interest. 

For  a  sample  calculation  of  attenuation  by  rain,  f  =  10  GHz  and  T  =  20  °C, 
figure  21  gives  the  result  that  cx  =  5  dB/km  when  R  =  150  mm/h.  At  a  range  of 
1  km,  az  =  5  dB  and  equation  (5)  gives  I(z)/I(z=0)  =  10-°-5  =  0.316.  At  ranges  of 
1  km,  an  intense  rain  is  seen  to  cause  attenuation  of  about  68  percent,  still 
less  than  an  order  of  magnitude. 


45 


Figure  16.  Attenuation  coefficient  a  of  rain  versus  frequency  for  f  =  0.1  —  10  GHz; 
compares  numerical  integration  of  the  MP  distribution  (T  =  40  °C)  with  discrete 
points  given  by  the  ORH  equation  (T  =  20  °C). 


46 


Figure  17.  Attenuation  coefficient  a  of  rain  versus  frequency  for  f  =  0.1  —  10  GHz; 
compares  numerical  integration  of  the  MP  distribution  (T  =  60  °C)  with  discrete 
points  given  by  the  ORH  equation  (T  =  20  °C). 
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8.  Summary 


In  this  report,  the  equations  which  describe  the  attenuation  of  long- 
wavelength  electromagnetic  radiation  by  rain  and  fog  were  briefly  reviewed 
and  explained  under  the  assumption  that  the  aerosol  droplets  were  spherical 
and  that  the  irradiance  of  the  beam  was  in  the  linear  regime. 

Calculations  of  water  fog  attenuation  coefficients  were  made  using  the 
Rayleigh  (i.e.,  small  particle)  approximation;  an  error  analysis  of  this 
approximation  was  made  by  comparing  Rayleigh  approximation  calculations 
of  absorption  efficiency  with  exact  Mie  theory  calculations  of  extinction 
efficiency.  Rayleigh  approximation  is  appropriate  for  fog  with  drop  diameters 
up  to  100  |iim  for  frequencies  less  than  50  GHz.  A  comparison  of  Rayleigh 
absorption  efficiencies  with  Mie  extinction  efficiencies  for  rain  with  water 
drops  up  to  5.0  mm  was  used  to  show  that,  for  frequencies  above  1  GHz,  the 
use  of  the  Rayleigh  approximation  for  rain  is  not  appropriate.  Numerical 
integration  of  the  Mie  extinction  efficiency  was  used  to  compute  the 
attenuation  coefficient  for  rain  with  an  MP  drop  size  distribution  for  various 
rain  rates  and  temperatures.  These  results  were  compared  with  those  given 
by  a  power  law  relationship  with  coefficients  given  by  ORH;  this  was  used  in 
the  NMMW  Module  of  the  EOSAEL.  The  ORH  relationship  gives  very  good 
results  at  temperatures  where  their  coefficients  were  determined  (from 
-20  to  20  °C).  As  temperature  increases  above  20  °C,  the  exact  Mie 
calculations  deviate  increasingly  from  the  results  of  ORH.  Above  20  °C,  the 
use  of  the  ORH  equation  is  not  recommended.  The  exact  Mie  computations 
may  easily  be  performed  for  any  temperature;  since  this  can  be  done  easily  on 
a  personal  computer,  this  approach  is  recommended. 
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Appendix  A 


Calculation  of  the  Complex  Index  of  Refraction 
of  Water  for  Frequencies  /  =  0  to  1000  GHz 
with  the  Double-Debye  Formulation 


The  calculations  of  attenuation  coefficient  a  done  in  this  report  as  a  function 

A 

of  frequency  require  that  the  complex  index  of  refraction  N  be  computed  as  a 
function  of  frequency.  In  the  body  of  this  report,  N  is  represented  for  a 
harmonic  wave  with  time-dependence  exp(-icot)  as  follows: 


N  =  Re(N)  +  ilm(N). 


(A-l) 


Liebe  et  al.  [3]  have  given  expressions  for  the  complex  permittivity  of  liquid 
water  e  (relative  to  the  air)  at  frequencies  f  =  0  to  1000  GHz  using  a 
double-Debye  formulation.  For  time-dependence  exp(-icot),  the  relative 
complex  permittivity  is  represented  as  follows: 

e  =  e'  +  ie.  (A-2) 

The  standard  relation  is  as  follows: 

N  =  (H£)’,  (A-3) 


where  the  magnetic  permeability  (J,  =  1  connects  the  real  and  imaginary 
components  of  N  with  those  of  e  and  the  result  obtained  is  as  follows: 


Re(N)  = 


e'  +  [(e')2+  (e")2] 


(A-4) 


hn(N)  = 


2Re(N)  ’ 


(A-5) 


The  absorption  efficiency  in  the  Rayleigh  approximation  is  noted  in  the  text  of 
this  report  to  be  proportional  to  Im(Z)  where  the  following  is  true: 


For  |i  =  1,  this  may  be  written  in  terms  of  the  relative  complex  permittivity  as 
follows: 


Z  = 


£-1 
£  +  2 


(A-7) 


For  calculations  with  a  hand  electronic  calculator,  it  is  useful  to  note  the 
following: 


Im(Z)  = 


3£a 


(e'  +  2)2  +  (£")2 


(A-8) 


The  double-Debye  formulation  of  Liebe  et  al.  [3]  gives  expressions  for  £'  and 
£a/  which  can  be  used  in  equations  A-4  and  A-7 : 


£'  -  (£0  -  £,) 


1  + 


\frJ 


+  (£]  -  £2) 


1  + 


Vfs7 


+  e2 


(A-9) 


£"  =  (£0  -  £j) 


1+ 


vfpy 


vfp  / 


+  (£1  -  £2) 


1+ 


ffV 

Vfs  7 


Vfs  7 


(A- 10) 


where 


£0  =  77.66+  103.3(6-1) 
300 


0  = 


£1 

£2 


°C 

=  5.48 
=  3.51 


-  +  273.15 


fp  =  20.09  -  142.4  (6  -  1)  +  294  (0  -  l)2  GHz 
fs  =  590 -1500  (0-1)  GHz. 
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Computations  of  Re(N)  andlm(N)  were  done  as  a  function  of  A,  for  T  =  0,  20, 
and  40  °C.  The  results  are  shown  in  figures  A-l  and  A-2.  Such  calculations 
were  earlier  performed  by  Falcone  et  al.  [6]  in  their  figure  13.  Comparison  of 
these  two  sets  of  results  shows  agreement  for  Re(N)  and  qualitative 
agreement  for  Im(N),  but  quantitative  differences  exist  for  Im(N).  This  is 
believed  to  be  caused  by  Falcone  et  al.'s  use  of  both  log  and  linear  scales  in 
their  figure  13.  [6] 


X  /  /xm 

Figure  A-l.  Real  part  of  the  complex  index  of  refraction  of  water  versus  wavelength 
for  X  =  102  —  10s  fim;  computed  using  the  double-Debye  expansion  with 
T  =  0  —  40  °C. 
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Figure  A-2.  Imaginary  part  of  the  complex  index  of  refraction  of  water  versus 
wavelength  for  X  =  102  —  10s  jim;  computed  using  the  double-Debye  expansion  with 
T  =  0  —  40  °C. 
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